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A New Catalog of Variable Stars in the Field of the Open Cluster 

M37 
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ABSTRACT 

We present a comprehensive re-analysis of stellar photometric variability in the field of the 
open cluster M37 following the application of a new photometry and de-trending method to 
MMT/Megacam image archive. This new analysis allows a rare opportunity to explore photo¬ 
metric variability over a broad range of time-scales, from minutes to a month. The intent of this 
work is to examine the entire sample of over 30,000 objects for periodic, aperiodic, and sporadic 
behaviors in their light curves. We show a modihed version of the fast periodogram algorithm 
(Fy^) and change-point analysis (CPA) as tools for detecting and assessing the significance of 
periodic and non-periodic variations. The benefits of our new photometry and analysis methods 
are evident. A total of 2306 stars exhibit convincing variations that are induced by flares, pulsa¬ 
tions, eclipses, starspots, and unknown causes in some cases. This represents a 60% increase in 
the number of variables known in this field. Moreover, 30 of the previously identified variables are 
found to be false positives resulting from time-dependent systematic effects. New catalog includes 
61 eclipsing binary systems, 92 multiperiodic variable stars, 132 aperiodic variables, and 436 flare 
stars, as well as several hundreds of rotating variables. Based on extended and improved catalog 
of variables, we investigate the basic properties (e.g., period, amplitude, type) of all variables. 
The catalog can be accessed through the web interface (http://stardb.yonsei.ac.kr/). 

Subject headings: binaries: eclipsing — catalogs — open clusters and associations: individual (M37) — 
stars: flare — stars: oscillations — stars: variables: general 


1. Introduction 


The rich open cluster M37 (NGC 2099) in 
the constellation Auriga has been photometrically 
monitored to search for new variables with dif- 


ferent motivations ( 

Kiss et al. 2001: Kane: et al.l 

2007; Messina et al. 

2008: Hartman et al.l 2008all. 


The combination of short and deep exposures with 
similar sky coverage allowed time-series photom¬ 
etry of hundreds to thousands of stars from the 
brightest to the dimmest within the cluster field 
(10 < V < 23). Their studies indicate that the 
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fraction of variable sources /var increases as pho¬ 
tometric precision of the survey gets better. 

Kiss et al. ( 200 ill performed the first variabil¬ 


ity survey for the central (29' x 18') field of the 
M37 with a i/c-filter photometry and found 7 vari¬ 
able stars among 2300 stars (/var ~ 0.3%). Due 
to poor photometric quality, their variables were 
with large variability and included three W UMa- 
type eclipsing binaries (EBs), two high-amplitude 
pulsating stars, and two long-period detached EB 
candidates. 

(l200Sh 


Kang et al.l (j2007h and iMessina et al 


used the same data set to search for short-period 
pulsating variables (e.g., S Scuti- and 7 Doradus- 
type stars) and rotating variables of later spectral 
types (F-K), respectively. Using U-band time- 
series observations, these two studies found the 
total of 153 variable stars among 12,000 stars 
(/var ~ 1.3%), including seven previously known 
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variables, in the field of the cluster (22.2' x 22.2'). 
Lastly, the most extensive var iability studies 


i n this field w ere carried out by iHartman et al 


( 2008bl . 2009al f9l as the by-products of a deep sur¬ 
vey for transiting planets in the M37 field (24' x 
24') with the 6.5m MMT telescope. The number 
of known variable stars increased by a factor of 
^10. By applying several statistical algorithms, 
they identified 1445 variables among 23,000 stars 
(■fva.r ^ 6.3%), of which 99% were new discoveries 
( Hartman et al. 2008bll . This catalog of variable 
stars provides the largest homogeneous set of stel¬ 
lar rotation periods for a cluster older than 500 
Myr, which showed the crucial role of the rota- 
tion in the stellar ang ular momentum evolution 
( Hartman et al. 2009a). 

In this paper, we search new v ariable phenom¬ 


ena from the same data set of Hartman et al 


( 2008afl with vastly improved photometric method 
( Chang et al. 20151 hereafter Paper I). In Section 
2 , we give a brief description of new light curves 
and its temporal properties. We present a com¬ 
prehensive re-analysis of photometric variability 
in Section 3. In Section 4, we present a new cata¬ 
log of variable stars and various variability charac¬ 
teristics. We compare our findings with previous 
works and discuss the nature of each variability 
class separately in Section 5. 


2. Data description 

We refer the reader to Paper I for a detailed 
description of our photometric reduction and light 
curve production applied to the archival imaging 
data of MMT transit survey. Here we provide only 
a brief summary of the key features. In Paper 
I, we introduced new methods for robust high- 
precision photometry from well-sampled images of 
a non-crowded held with strongly varying PSF. 
The archival light curves from the original im¬ 
age subtraction procedure exhibit many unusual 
outliers, and more than 20% of data get rejected 
by the simple hltering algorithm adopted by early 
analysis. We note that some of these outliers are 
results of intrinsic variability of short time scales. 
In order to achieve better photometric precisions 
and also to utilize all available data, the entire 
imaging database were re-analyzed with our multi¬ 
aperture indexing photometry and a set of sophis¬ 
ticated calibration procedures. This algorithm 
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Fig. 1.— Observational cadence of the MMT tran¬ 
sit survey, showing cumulative fraction function 
of time intervals ATmjd between two consecutive 
frames (top left) and the exposure times Texp (top 
right). We use logarithmic scale for the horizontal 
axis in these two panels. The bottom panel is the 
distribution of exposure times as a function of the 
modified Julian date (MJD). 


finds an optimal aperture for each star with a max¬ 
imum signal-to-noise ratio, and also treat pecu¬ 
liar situations where photometry returns mislead¬ 
ing information with more optimal photometric in¬ 
dex. We also applied photometric de-trending al¬ 
gorithm based on a hierarchical clustering method. 
Our method removes systematic variations that 
are shared by light curves of nearby stars while 
true variabilities are preserved. Consequently, our 
method utilizes nearly 100% of available data and 
reduce the RMS scatter several times smaller than 
original analysis for brighter stars. 


Our new data set includes 30,294 objects, in¬ 
cluding both point and extended sources. For 
blended targets, the original image subtraction 
light curve is still of higher precision than the 
present aperture photometry light curves, and 
thus for these objects the archival light curve is 
used. In order to pick out such cases, we utilized 
a simple method for quan tifying the level of blend¬ 
ing (see Irwin et M]l2007h . 
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where 


(1) 

p 

for light curve points nip with uncertainties Up, 
and Xpoiy i® the same statistic measured with re¬ 
spect to a 4th order polynomial fit. We adopt the 
value b > 0.6 for selecting significant blended light 
curves, and then count the number of blended ones 
(3746 objects). 

Figure [T] shows that the cumulative fraction 
function of temporal intervals AT^jd between two 
consecutive frames and the exposure times Tg^p- 
This observing cadence is sensitive to variability 
on timescales ranging from tens of seconds (30 
seconds) to one month (^31.7 days). In partic¬ 
ular, the time resolution allows an analysis of the 
short-term variability, which is less than ~80 sec¬ 
onds over 50% of the data and ~150 seconds over 
90% of the data. 


3. Variability-analysis methods 


Variability features can be divided into two 
broad categories: periodic and non-periodic fea¬ 
ture. The former is easier to detect and model 
using statistical and model-specific fitting proce¬ 
dures as a form of Fourier series (e .g. JPebosscher et al 
2007 . 20091 : iRichards et al.l 2011 1. But the latter 
is rather ambiguous, since it is not always possible 
to characterize flux variation and to define a reli¬ 
able rang e for some parameters. To overcome this 
difficulty, [Richards et al.l l|201l[) calculated many 
basic statistics that can describe the distribution 
of fluxes even in the limit of few data points. 
One interesting approach is to select stochasti¬ 
cally varying features in quasar-like sources. The 
correlated variability of quasars is well parameter¬ 
ized by a damped random walk model, providing 
a simple statistical description of thei r character¬ 
istic timescales and amplitu des (e.g., iKellv et "ni 
20091 : Butler fc BloomI 20 111) . In this work we use 
a modified version of the fast periodogram algo¬ 

rithm (Fy^) and the change-point analysis (CPA) 
method for detecting and assessing the significance 
of periodic and non-periodic signals, respectively. 


3.1. Generic periodic feature detection 

We first searched a periodic feature by Fy^ 
technique that calculates the minimized y^ as a 
function of frequency at the desired number of 


harmonics ( Palmer 20091) . With the best-fit fre¬ 
quency, we compared the ratio between the vari¬ 
ance of the data after subtracting the best-fit 
model and full-amplitude of model function. This 
algorithm use a Fourier series truncated at har¬ 
monic h as a model function: 


= Ao -I- ^ A 2 h-i sin(/i27r/t) -|- A 2 h cos(/i27r/t) 
h=l 

( 2 ) 

The fitting coefficients (Aq, A 2 / 1 , and A 2 ?i_i) ob¬ 
tained in each harmonic model provide a reason¬ 
able description of a generic periodic feature, ir- 
respective of the featur e shape, as discussed by 
Debosscher et al. ( 20091) . Figure [5] shows exam¬ 
ple light curves of both periodic (new variable 
star; V1983) and long-term (V1092) variable stars 
in our sample, respectively. Because the optimal 
choice of h depends on the true shape of the sig¬ 
nal, the fit quality is evaluated for each model 
with standard y^ statistics, and then checked the 
residual (phase-folded) light curves after the sub¬ 
traction of the best-fit model. To test the sig¬ 
nificance of a signal, we use the ratio between 
the full-amplitude of the model function (Famp 
= nT-max —mmin) and the RMS scatter {aa) of the 
data after model subtraction: 


Famp > Xaa- (3) 

The lower limit for the ratio is empirically set to A 
= 0.95 based on light curve inspections. A-values 
less than 0.95 indicate that we do not expect any 
significant (semi-) periodic signals for a given light 
curve. Even when the source is poorly described 
by the harmonic fit, this indicator provides a hint 
of a generic periodic signal (see right panel of Fig¬ 
ure [2]). 

3.2. Arbitrary-shaped and sporadic event 
detection 

3.2.1. Change-point analysis (CPA) 

Most variability detection methods require con¬ 
ventional models that are mainly focused on the 
periodic signals, and are not suitable for the study 
of arbitrary-shaped, non-periodic, and sporadic 
variations. Also, in many cases, signal estima¬ 
tion is equated with smoothing of data for de- 
noising. This sometimes discards vital informa¬ 
tion in time series data. We introduce a non- 


3 

































0 0.25 0.5 0.75 1 0 0.25 0.5 0.75 

Phase 


0.25 0.5 0.75 1 0 0.25 0.5 0.75 1 

Phase 


Fig. 2.— Example light curves of a periodic star (four plots at left panel) and another quasi-periodic star 
(four plots at right panel) with different periods detected by the first feature detection algorithm, respectively. 
The shape of periodic feature is varied widely by the selected harmonic model h (from 1 to 4), but we can 
choose the best candidate with statistics. 


parametric method to extract all significant fea¬ 
tures based on the CPA using local statistics. 


Change-point analysis is a method to identify 
abrupt variations in sequential data. It is widely 
used in the statistics and data mining communities 
as well as in the field of time domain astronomy 
(see e.g., LSchiitz fc Holschneidei 2011 and refer¬ 
ence therein). Using a combination of cumulative 
sum scheme (CU SUM) and bootstrap rank statis¬ 
tics (lTavloi<20nnll . we produce a series of estimated 
change-points which correspond to the moments of 
apparent systematic changes. For a given dataset 
{xi, X 2 , ■ • ■, Xn} of size n, the CUSUM values are 
given by 


t 

St = y^(xi - x) 


(4) 


i=l 

for t = 0,1, • • • , n, where 5'o = 0 and the mean 
value X for the sample of n values of Xi. An in¬ 
flection point at which the sign of the CUSUM 
slope changes is used to determine whether a 
given interval of data should be kept as one 
{xi=X2=‘ • • =Xn=x) or subdivided into two subin¬ 
tervals (xi=- ■ ■ =Xp=/^Xp+i=Xn)- Based on N{> 
1 ,000) bootstrap samples that are randomly re¬ 
ordered original values, we estimate the confidence 
level (c.l.) to reduce the false positives by random 
noise. Here are the main procedure of each boot¬ 


strap process: 

• We generate a bootstrap sample of n dataset, 
denoted {x\, x\, • • •, a;^}, by randomly re¬ 
ordering the original n values without re¬ 
placement. 

• Based on the bootstrap dataset, we calculate 
the bootstrap CUSUM by using equa¬ 
tion (4), and then obtain the maximum, 
minimum, and difference of the bootstrap 
CUSUM. 


qb — qb _ qb 
^diff ^max 


(5) 


where 


5(^,, = max{.So^5^•■• ,50}, 

5L„=min{50,5},.-- ,50). 


• We check whether the bootstrap difference 
Sdiff i® 1®®® than the original difference 
Sdiff — Smax Sjjiiji . 


It is straightforward to derive estimates of confi¬ 
dence level as follows: 

c.l. = 100 X (6) 
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Fig. 3.— Example light curves of flare-like (left panel) and eclipsing-like (right panel) events detected by the 
second feature detection algorithm. We used the default values of the variable parameters (W=3, iV 2 =l, 
iV3=3) as the detection thresholds. The CPA model is over-plotted in dashed lines. 


where 




1 if ^diff < ^diff 

0 if ^dzff ^ ^diff 


characterized by the start and end time of a given 
interval with its mean and variance. After this 
procedure, signal detection process is performed 
for each segment separately. 


Typically more than 90% conhdence level is suit¬ 
able for most cases. When a sudden change is 
detected, the location of change-point is initially 
determined as follow: 

Pt = argmax|5't|, (7) 

iG[0,n] 

where pt denotes the last point before the change 
occurred. The time-series data is split into two 
segments on each side of the change point, and 
the analysis is repeated for each segment until no 
more significant change point is detected. 

Once this set is generated, all the change- 
point locations and their confidence levels are re- 
estimated in a backward elimination manner. Us¬ 
ing only the sample in the segment bounded by 
adjacent change-points, we re-calculate these pa¬ 
rameters detected by the forward procedure to ob¬ 
tain more accurate estimations. When any initial 
change-point fails the criteria in the backward pro¬ 
cedure, we eliminate this change-point estimate. 
This process repeats until no further improvement 
is possible. The final change points thus define 
the segments (i.e., piecewise constant level sets). 


3.2.2. Example of flare-like or eelipse-like event 
detection 


To detect significant outlying features occur¬ 
ring at specific levels in light curve, we adopt a 
simp le criteria similar to Micro-lensing Alert sys¬ 
tem (|Glicensteiril 1200 ill in the presence of hetero- 
scedastic measurement errors {wi)\ 


Xi- XL < 0, 

(8a) 


(8b) 

CTL 

\Xi -XL+Wi\ 

> t^2, 

(8c) 

UL 

ConM > Ns, 

(8d) 


where the mean xl and deviation are the lo¬ 
cal statistics for a given segment, Wi is the pho¬ 
tometric error at epoch i, and ConM is the num¬ 
ber of consecutive points which satisfy the equa¬ 
tions (8a-8c). The different combination of Afi, 2,3 
values are used for maximize the detection effi¬ 
ciency of significant outlying features. The values 
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of -/Vi. 2,3 were taken to be at least larger than 3, 
1 , and 3, respectively. As a result, statistically 
significant events were flagged as candidate flare 
sources if at least three contiguous points satisfy 
the above criteria. Example of flare-like feature 
detection is shown in the left panel of Figure [3] 
By applying the CPA method to the entire light 
curves, we efficiently identified several hundred in¬ 
stances of abrupt brightness changes without any 
smoothing or interpolation of the raw data. 

We use the same limits of CPA thresholds to 
find eclipsing-like event detection, except the first 
term of equation (8), i.e., Xi — xl > 0. The 
CPA method is useful to detect the moments of 
eclipse ingress, center, or egress in cases where 
the eclipsing pattern is not repeated or the cover¬ 
age is not sufficient for the detection through con¬ 
ventional period analysis (e.g., Box-fitting Least- 
Squares method). Example of eclipsing-like fea¬ 
ture detection is shown in the right panel of Fig¬ 
ure [3l We will discuss this issue further in Section 
5.1. 


3.3. Periodicity analysis 

For periodicity, we first used the VARTOOLS 
light curve analysis program to determine likely 
periods of variable star|3. This code com¬ 
putes three different period ograms for the light 
curves;(i) Lomb-Scargle (LS: iLomb 1976 Scargld 


Il982l:l^ess fc Rvbickill98^|Press et al.lll992lf 

^oV:E 


Anal y sis of Varianc e ( AoV: ISchwarzenberg-Czernv 
[iiii [5 evo 2 l2nn,5h. and fiii) B ox-fitting Least- 
Squares lBLS: lKovacs et al ■I2nn2lf algorithms. Be¬ 
cause each method is sensitive to different types 
of periodic variables, it is often useful to combine 
results from many different analysis of periodic- 
it y. Basically we adopt t he similar treatment as 
in iHartman et all ( 2008bh . We limited the search 
range of these algorithms to periods between 0.01 
and 30 days with relatively fine frequency resolu¬ 
tion of less than 0.01/T (T is the time-span of the 
light curve). 

We also obtained period estimates from the sta¬ 
tistical technique SigSpec ( Reegenir2007 \ It cal¬ 
culates the spectral significance (sig) levels based 
on the analytical solution for discrete Fourier 
transform amplitude spectrum, including depen¬ 


^http: //www. astro .princeton. edu/~ jhartman/vartools .html 
^http://homepage.univle.ac.at/peter.reegen/download.html 


dencies on frequency and phase. In recent years, 
this package is mainly used in the study of stel¬ 
lar pulsation (i.e., asteroseismology), which shows 
up as small, periodic variations in the bright¬ 
ness of the stars fe.g.. [&uberbauer et al.j2007 : 
Benko et al. I l201(ll: iGuggenberger et ~ I l2nl2^ . In 


order to further crosscheck results of periodic¬ 
ity analysis, we also performed a multi-step pe¬ 
riod search algorithm MS_Period ( Shin fc Bvun 

2on4li^ . 


Figure ID shows examples of three selected light 
curves; power spectrum and phased diagram cor¬ 
responding to the highest peak. The differences 
in the various algorithms depend strongly on the 
shape of variation. One case (V2256) shows that 
all periods are essentially identical {P ~ 0.0784 
days), while in the other cases the derived period 
differ a lot. Visual inspections are necessary for 
V2173 and V1482 to choose the most plausible 
periods; P ~ 0.726 and 2.595 respectively. 


4. New catalog of variable stars in M37 
field 


We provide a catalog of the total 2306 stars 
that exhibits convincing variations that are in¬ 
duced by flares, pulsations, eclipses, starspots, 
or in some cases, unknown causes. This cat¬ 
alog is made available in electronic form at 
http: //stardb. yonsei. ac. kr/ Table lAfl gives 
the basic information about previously known and 
new variable stars that we detected in the field of 
M37. The columns are as follows: 


VarlD: The variable star IDs for the first 


1483 stars are from the numbering system of 
Hartman et al. ( 2008bl l. New variable ID num¬ 


bers are given in ascending order by their right 
ascension from V1484 to V2306. 


Object designation: Each object is designated 
by the combination of right ascension and declina¬ 
tion for J2000 epoch (e.g., 055536.78-k323345.97). 

Period : We basically used the dominant peri¬ 
ods which were highly scored by several methods. 
For EBs, we checked the existence of secondary 
eclipse when phased at twice the period. In cases 
where our periodicity analysis is affected by a large 
scatter of light curves due to photometric bias or 
source blending, we used periods taken from the 


^ http: //www. physics. ox. ac.uk/users/nisshin/scieiice/code/MultiStep_Period/ 
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Fig. 4.— Examples of power spectrum measurement and its corresponding phased diagram with period at 
the highest peak (arrows) in power spectrum. The scales of y-axis are normalized to the maximum value. 
One case shows that all periods are identical within errors (left panel; V2256), while in the other cases the 
derived period differs a lot. 


previous catalog ( Hartman et al.l[2008b ). and then 
denoted by f hr that column. We also found evi¬ 
dences of multiperiodicity in some pulsating stars 
(see Section 5.2). 

Amplitude: The range of variability for peri¬ 
odic variables was derived from phase diagram 
in which the shape of the light curve was mod¬ 
eled by equation (2). This amplitude is defined 
as the full-amplitude of model function (Amp = 
nT-max — ’mmin)- When this was not possible, i.e., 
cases of EBs and irregulars, we used the original 
light curve to estimate the amplitude. 


Magnitudes: ugrizJHKg magnitudes for many 
of our stars are available from the Sloan Digital 
Sky Survey Data Release 7 (SDSS DR7jl and Two 
Micro All Sky Survey (2MASS)11. Each catalog is 
positionally matched to MMT sources with an av¬ 
erage distance of 0.6". In the case of SDSS, the 
photometric quality at low Galactic latitudes is 
not guaranteed to be accurate to the SDSS quoted 
limits of 2% in color, and photometric depths are 
about two magnitude shallower than our data. 
Similarly, 2MASS is not deep enough to detect 
low-mass members in the cluster. 


VarType: We assigned five variability types 
similar to the variability classification scheme pro- 


^http://www.sdss.org/dr7/ 

®http://www.ipac.caltech.edu/2mass/ 


posed by the GCVS ( Samus et al.l2009ll . Variable 
stars are grouped according to major astrophysical 
reasons for variability, such as flare (F), pulsating 
(P), rotating (R), eclipsing binary (EB), and ape¬ 
riodic variables (A), including variable candidates 
(:). Further details are given in Section 5. If a vari¬ 
able belongs to more than one type of variability, 
it is indicated by combination of individual types, 
e.g., F+EB, P+EB:. 

A total of 30 sta rs listed as variables in 
Hartman et al. ( 2008bl) turned out to be not real 


variable in our study. Their new amplitude spec¬ 
trum often does not show anything signihcant 
(Figure[5]) and no apparent periodicity in the light 
curve. These false positives are listed in Table [H 


4.1. Variability fractions 

In this work, we found 823 new variables 
(V1484“V2306) and identified types of their vari¬ 
ability. The discovery rate of new variables is 
increased by 60% in comparison with the exist¬ 
ing catalog. The new variables are primarily 436 
flaring stars, 30 sharp EBs, 48 low-amplitude pul¬ 
sating variables, and 65 aperiodic variables, as 
well as many other rotating variables. 

We briefly explored how variability fractions 
depends on their brightness/period and variable 
types. As shown in Figure [51 short-period, low- 
amplitude pulsating stars and episodic flare events 
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Table 1 

List of 30 stars whose variability turn out to be false 


No. 

VarlD 

StarlD 

(J2000) 

Chip 

(#) 

r 

(mag) 

Period^ 

(days) 

1 

V27 

055118.34+323201.71 

32 

20.938 

16.364 

2 

V85 

055125.32+324341.29 

28 

20.862 

7.668 

3 

V88 

055125.46+324315.32 

28 

21.820 

1.646 

4 

V120 

055128.52+322339.10 

35 

17.196 

6.859 

5 

V155 

055131.10+323756.98 

30 

16.847 

7.424 

6 

V166 

055131.96+323600.24 

31 

20.895 

10.451 

7 

V224 

055137.29+324011.78 

29 

18.165 

4.571 

8 

V230 

055138.02+324011.65 

29 

16.690 

4.595 

9 

V250 

055139.59+324013.27 

29 

17.282 

4.551 

10 

V330 

055151.25+323527.70 

22 

16.831 

5.428 

11 

V337 

055151.72+323510.02 

22 

17.920 

5.442 

12 

V342 

055152.18+323522.44 

22 

19.061 

5.438 

13 

V347 

055152.55+323536.03 

22 

15.574 

5.500 

14 

V367 

055153.88+323459.77 

22 

17.463 

5.515 

15 

V390 

055155.42+323523.56 

22 

18.268 

5.481 

16 

V408 

055156.77+323534.61 

22 

15.733 

5.481 

17 

V433 

055157.99+323528.59 

22 

16.617 

5.481 

18 

V492 

055202.16+324517.30 

19 

20.849 

16.586 

19 

V973 

055235.48+322104.58 

18 

16.847 

5.089 

20 

V1076 

055241.54+322107.71 

18 

20.262 

10.927 

21 

V1237 

055258.14+322152.09 

9 

16.775 

5.692 

22 

V1260 

055300.53+323445.47 

4 

19.851 

15.895 

23 

V1282 

055303.05+323621.95 

4 

21.530 

1.666 

24 

V1305 

055305.05+323447.83 

4 

18.166 

7.966 

25 

V1330 

055306.74+323442.52 

4 

18.061 

8.143 

26 

V1371 

055310.61+323748.13 

3 

20.902 

9.934 

27 

V1372 

055310.65+323740.82 

3 

20.985 

10.633 

28 

V1400 

055312.89+323005.46 

6 

16.770 


29 

V1421 

055314.64+322656.12 

7 

20.735 

15.013 

30 

V1423 

055315.29+323051.00 

6 

20.728 

12.388 


Period is obtained bv IHartman et all Il2008bfi . 












Fig. 5.— Amplitude spectrum of the archival 
(gray) and our new light curves (black). Mean 
noise levels are greatly reduced in new analysis. 


are dominantly found at the bright- and faint- 
end of magnitude, respectively. In Figure [71 we 
plot variability fraction as a function of period. 
The period range of new variables is very si milar 
to that reported by [Hartman et ah ( 2008bll . Of 
these, most of new pulsating candidates show 6 
Scuti-like characteristics {P < 0.3 days). Also the 
majority of periodic variability can be explained 
by rotational modulation of surface features (i.e., 
starspots). 


4.2. Time-scales and amplitudes of the 
variability 

Time-scales and amplitudes of variability are 
intrinsic proper ties of the star and are independent 
of the distance ( Ever fc Mowlavill2008ll . These two 
quantities can easily be measured from the light 
curves. We extracted the basic properties (e.g., 
period, amplitude) of all periodic variables using 
the Fourier decomposition technique. For compar¬ 
ison purposes, we utilized the flare parameters: (i) 
flare duration to. 2 (= to, 2 ,nse+to. 2 ,decay) that is de¬ 
fined as the time-scales of flare events at 80% of its 
peak level, and (ii) flare amplitude AiUpeak that is 
magnitude at the moment of flare maximum (S. 
-W. Chang et ah, in preparation, hereafter Paper 
III). 

We separate the variable objects into photomet¬ 



10 12 14 16 18 20 22 24 26 

r(mag) 


Fig. 6.— Variability fraction as a function of 
r magnitude. Previously known variables (VI- 
V1483), newly detected variables (V1484-V2306), 
and total number of observed stars are marked 
with dashed lines, gray solid lines, and black solid 
lines, respectively. 


rically selected c andidate cluster and fie ld mem¬ 
bers provided by Hartman et al.l ( 2008bll . respec¬ 
tively. Since the g — i color covers a relatively wide 
range from —0.5 to 4.5, we used g — i as a. sepa¬ 
rator between a sequence of cluster and that of 
field population. As can be seen in Figure [H short 
period (< 1 day), small amplitude pulsating stars 
are found at the blue-end of the distribution of 
both cluster and field sources, while episodic flare 
events are dominantly found from the sources in 
the red color range. 


Our new data strengthen the relationship be¬ 
tween the color of stars and their variability pa¬ 
rameters derived for cluster members previously 
(Figure [8]). The cluster members show a tight 
and nearly linear relation between color index and 
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Fig. 7.— Variability fraction as a function of pe¬ 
riod. The representation of the lines is same as 
Figure m 


rotation period in the range of 1 < logP < 20 
days, although the change in slope occurs after the 
knee point at P ~ 7 days. This sequence forms a 
curved, diagonal band from bright blue stars to 
faint red ones, whose periods are increasing with 
increasing color index (or magnitude), while for 
non-cluster members there is no clear correlation. 


The color-amplitude relation for the cluster 
members shows that, statistically, cooler stars 
have greater amplitudes in the case of both pe¬ 
riodic variables and flares. As in the previous 
study on stellar rot ation (see Section 6 an d Fig¬ 
ures 14, 15 and 17 in Hartman et al.l[^09all . there 
is a strong correlation between Rossby number and 
amplitude, which can be interpreted as stars rotat¬ 


Fig. 8.— g — i color vs. variability time-scales 
and amplitude relations for all periodic variables 
(black points) in the cluster and field region, re¬ 
spectively. For comparison purposes, we also show 
the parameter distributions of 625 flare events 
(gray squares). The flare time-scale is defined as 
the time duration of flare events at 80% of its peak 
level (T 0 . 2 ), and the flare amplitude is magnitude 
at the moment of flare maximum. 


ing faster having more spots, larger spots, and/or 
greater spot-photosphere contrast than slower ro¬ 
tating stars. Cooler stars have on average lower 
Rossby numbers at a particular age, so they end 
up having larger amplitudes. 

4.3. Color-magnitude diagram of variable 
stars 

Figure |9] illustrates the color-magnitude dia¬ 
gram (CMD) for all objects and all variable stars 
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5. Types of variability 



24 


-1 012345 -1 012345 

g-i 9-i 

Fig. 9.— Color-magnitude diagram with marked 
positions of variables (black dots) found in the 
M37 held. We compare the g — i vs. i diagram 
of both observed (left) and simulated (right) data. 
Photometrically selected cluster region is plotted 
as red line in both panels. 


in the M37 Held. In the new catalog, we found 
that about half of variable objects are distributed 
along the cluster main-sequence, and the other 
half are located in the Held or in the outskirts 
of the cluster sequence. In order to identify the 
latter stars, we compared our data with the TRI- 
LEGAL@ model to reproduce the CMD at this 
location of the Galaxy. Most of the input pa rame- 

(I 2 OO 5 II . For 


ters are described in the Girardi et al. 


this analysis, we set the pointing parameters for 
the Held center (Z = 177.635; b = 3.091) and total 
Held area (□ = 0.16 deg^), a limiting magnitude as 
r = 23.5 and a Salpeter initial mass function with 
cutoff at 0.01 M 0 . The remaining parameter val¬ 
ues are fixed to default values that correspond to 
their model calibration. The model output con¬ 
sists mostly of thin disk stars (96%) with small 
contribution of old stars in the thick disk and halo 
(4%). Non-member, variable stars may be related 
in some way to the thin disk population, but their 
distributions are not fully explained by the model. 


®The latest version of TRIdimensional 
modeL of thE GALaxy is available at 
http://stev.oapd.inaf.it/cgi-bin/trilegal 


In the following subsections, we investigate the 
nature of each variability class: (i) 61 EB systems, 
(ii) 92 multiperiodic variable stars, (Hi) 132 ape¬ 
riodic variables, and (iv) 436 Hare stars, respec¬ 
tively. In order to understand variability proper¬ 
ties, we utilized both color-magnitude and color- 
color diagrams to separate variable objects into 
broad magnitude ranges or spectral classes, to the 
extent possible. 

5.1. Eclipsing binary star systems 

Our photometric time-series data provide the 
light curves of 61 EBs. Figure [TO] shows phased 
light curves of 8 newly discovered EBs. We catego¬ 
rize EBs into three groups according to the shape 
of phased light curves, and discuss each group be¬ 
low. 

5.1.1. Detached or Semi-detached EBs 

EBs in this category are either Algol-type sys¬ 
tems or P Lyrae-type systems. To derive their ge¬ 
ometric parameters, we used the Detached Eclips¬ 
i ng B inary Light curve (DEBiL) fitter by IPevor 
(12005^ ^. which employs a simple model of a per¬ 
fectly detached binary system in a classical 2-body 
orbit. Table [2] lists parameters derived for all 
detached/semi-detached binaries in the M37 field. 
EBs with well determined periods are included in 
DEBiL analysis. Special notes are given below for 
interesting features. 

V2 This bright EB system was missed by pre¬ 
vious study due to magnitude cut, but it clearly 
shows evidences of eclipse in our new light curve. 

V587 The light curve of this EB system con¬ 
tains two distinct periodic signals, the 3.6173 days 
modulation with a small amplitude (Am ~ 0.056 
mag) and an additional period at 3.2857 days (or¬ 
bital period). Since an ellipsoidal variation should 
have the same period as the orbit, the first signal 
is more likely a consequence of starspots and ro¬ 
tational modulation.. After 3.6173-d modulation 
is removed, the light curve clearly shows the sig- 


^https://www.cfa.harvard.edu/-jdevor/DEBiL.html 
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Fig. 10.— Phased light curves of 8 newly discovered eclipsing binaries. 


nature of a typical Algol-type binary with both 
primary and secondary eclipses. 


V1482 This EB system is a newly recognized 
Algol-type binary with a 6 Set-type pulsating com¬ 
ponent. We have clearly detected an orbital pe¬ 
riod of 2.5972 days including multiple primary and 
secondary eclipses. The short period variability 
caused by stellar pulsation is shown in the out-of¬ 
eclipse light curves. In order to determine the os¬ 
cillation modes of the pulsation, a frequency anal¬ 
ysis was performed in the o ut-of-eclipse data usin g 
the program PERIOD04 ( Lenz fc Breged[2005 1. 
Our analysis not only shows peaks as a conse¬ 
quence of binary effects (/i = 2.5638 and /2 = 
1.3006) but also reveals the presence of the domi¬ 
nant pulsation frequency /a = 0.0465 with a small 
amplitude (Am = 0.0035 mag). 


5.1.2. Contact EBs 

The light curves of contact binary stars (also 
known as W UMa-type variable stars; EW) are 
very much different from those for detached bi¬ 
naries. Our new catalog lists 22 contact binaries 
of which four are new discoveries (V840, V1031, 
V1250, V1463). 


V37, V706, 


V1031 These systems are one of 
the fastest orbiting contact binary systems (P ~ 
0.22 days). It is known that observed distribu¬ 
tion of orbital periods has a cut-off p eriod limit at 

120021 [ 


Molnaretal 


0.22 days fe.g.. lRucinski 

I 2 OI 3 II . We have not found anything faster in our 
study. 


5.1.3. New EB candidates with eclipsing-like fea¬ 
tures 

With our CPA analysis, we discovered 15 new 
eclipsing binary candidates that have eclipsing-like 
or transit-like features undergoing eclipse, eclipse- 
ingress, or eclipse-egress phases. Table [3] summa¬ 
rize observed eclipsing features of the light curves, 
which are missed by the standard analysis soft¬ 
ware. Their orbital periods are not well con¬ 
strained. The light curve of each event can be 
described by the first time of dimming points to, 
the time duration between the start time of eclipse 
and the end time of eclipse At, and the depth of 
observed dimming points d. Figure [TT] illustrates 
examples of light curves of this category. 


V791 This detached EB system has deep pri¬ 
mary and secondary eclipses (12-13% and 1%, re- 
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Table 2 

Parameters for 22 Detached and Semi-detached EB systems 


No. 

VarlD 

r 

(mag) 

Porb 

(days) 


Ri/a 

R 2 j a 

mi 

(mag) 

m 2 

(mag) 

sin(2) 

+ b 

to 

UJ 

(deg) 

Note^ 

1 

VI 

13.677 










Not enough data 

2 

V2 

14.840 










Not enough data 

3 

V55 

22.549 

0.7780 

<0.08 

0.17 

0.15 

22.81 

24.12 

1.00 

0.65 

283.4 

D 

4 

V146 

20.463 

0.6477 









SD; No secondary minimum? 

5 

V226 

22.121 

0.9123 

<0.01 

0.10 

0.04 

22.23 

24.59 

1.00 

0.53 

129.5 

D 

6 

V395 

17.392 

3.4460 

0.01 

0.10 

0.08 

17.96 

18.35 

1.00 

0.33 

232.1 

D 

7 

V457 

19.474 

2.1995 

<0.08 

0.08 

0.07 

20.15 

20.30 

1.00 

0.75 

239.7 

D 

8 

V485 

17.080 

1.6676 









D; No secondary minimum? 

9 

V587 

16.572 

3.2857 

0.52 

0.06 

0.02 

16.58 

21.24 

1.00 

0.31 

209.4 

D -h RM (P=3.6173 days) (new) 

10 

V692 

16.669 

1.8324 

0.06 

0.28 

0.06 

16.73 

19.80 

0.98 

0.01 

77.2 

SD 

11 

V716 

21.635 

2.0272 









D; No secondary minimum? 

12 

V733 

17.383 

0.7310 

<0.04 

0.43 

0.34 

17.61 

18.91 

0.99 

0.77 

73.3 

D 

13 

V965 

16.736 

1.9732 

0.00 

0.16 

0.05 

16.75 

21.05 

1.00 

0.27 

357.0 

SD (new) 

14 

V1028 

18.730 

5.4956 

0.03 

0.04 

0.04 

19.21 

19.84 

1.00 

0.55 

267.5 

D + RM 

15 

V1141 

16.977 

3.6948 

< 0.66 

0.06 

0.01 

16.98 

22.60 

1.00 

0.29 

268.3 

D 

16 

V1182 

21.869 

2.2834 

0.05 

0.07 

0.06 

22.40 

22.87 

1.00 

0.16 

268.3 

D (new) 

17 

V1187 

19.742 

1.1919 

<0.41 

0.12 

0.02 

19.77 

23.56 

1.00 

0.86 

271.6 

D 

18 

V1380 

14.874 

2.1898 

<0.05 

0.16 

0.06 

15.01 

17.19 

0.99 

0.84 

280.1 

D -h 7 Doradus component (P=0.941/1.174 days) 

19 

V1456 

17.006 

5.9680 

0.12 

0.07 

0.07 

17.47 

18.15 

1.00 

0.54 

89.3 

D 

20 

V1482 

15.597 

2.5972 

0.02 

0.14 

0.08 

15.71 

18.14 

0.99 

0.17 

76.1 

D -h (5 Scuti component (P=0.0465 days) (new) 

21 

V2138 

14.181 

3.4051 

0.04 

0.09 

0.04 

14.28 

16.85 

1.00 

0.63 

191.7 

D (new) 

22 

V2297 

15.863 

3.5109 

<0.01 

0.12 

0.02 

15.88 

20.42 

1.00 

0.32 

229.7 

D (new) 


Note.— The column description is as follows, e: Orbital eccentricity; i?i/a: Radius of large star (in units of semimajor axis); R 2 / 0 : Radius of small star (in units of 
semimajor axis); mi: Brightness of large star (magnitudes); m 2 : Brightness of small star (magnitudes); sin( 2 ): Sine of inclination; to: Phased epoch of periastron; lj: 
Argument of periastron (in degrees). 

®^For the cases where the errors are larger than e, we show the best-fit value plus the error as an upper limit. 

^Phased epoch of periastron: Modified Julian date folded by the period (see Table 4 in |DevoHl2005D . 

detached eclipsing binary system; SD: semi-detached eclipsing binary system; RM: rotational modulation. 








Table 3 

List of new 15 eclipsing binary candidates from CPA analysis 


No. 

VarlD 

Star 

(J2000) 

r 

(mag) 

^event 
(#) 

io" 

(day) 

aA 

(hour) 

(mag) 

1 

V791 

055223.89+322826.47 

15.604 

1 

53732.28145 

1.0 

0.011 





2 

53734.26454 

2.6 

0.142 





3 

53756.09692 

3.1 

0.134 

2 

V827 

055226.24+324044.62 

15.774 

1 

53731.14064 

1.1 

0.074 





2 

53745.26403 

4.7 

0.171 

3 

V926 

055232.65+322208.16 

16.348 

1 

53755.08741 

2.2 

0.046 

4 

V985 

055236.13+323001.13 

16.012 

1 

53734.15888 

2.2 

0.054 





2 

53756.24427 

1.1 

0.051 

5 

V1018 

055237.96+323415.64 

16.569 

1 

53732.14053 

1.5 

0.019 





2 

53756.17928 

1.8 

< 0.030 

6 

V1483 

055311.98+322353.58 

17.717 

1 

53726.20141 

0.4 

0.208 





2 

53730.29695 

0.7 

0.146 





3 

53733.22359 

1.6 

0.184 





4 

53734.41413 

0.4 

0.159 





5 

53737.34213 

1.2 

0.161 





6 

53740.33536 

0.9 

0.124 





7 

53753.25941 

0.7 

0.203 





8 

53756.17828 

1.2 

0.157 

7 

V1731 

055157.49+324534.90 

15.823 

1 

53728.25931 

3.3 

0.233 





2 

53732.11914 

2.5 

0.282 

8 

VI748 

055159.96+322433.64 

14.846 

1 

53740.46821 

0.5 

0.012 





2 

53742.46911 

0.5 

0.015 

9 

V1757 

055201.06+322408.79 

17.189 

1 

53731.14064 

0.4 

0.022 

10 

V1889 

055216.11+322822.95 

14.971 

1 

53742.42521 

1.6 

0.100 





2 

53756.13038 

3.8 

0.082 

11 

V1892 

055216.39+322609.00 

19.328 

1 

53734.15888 

2.0 

0.178 

12 

V1909 

055218.62+324100.16 

13.281 

1 

53731.34461 

0.0 

0.019 





2 

53753.15826 

2.0 

0.291 

13 

V1915 

055218.84+324434.08 

21.189 

1 

53755.21194 

1.9 

0.273 

14 

V2213 

055301.96+323744.91 

16.433 

1 

53725.17392 

2.0 

0.230 





2 

53733.15900 

0.5 

0.025 

15 

V2214 

055302.02+322327.62 

15.094 

1 

53732.11914 

2.6 

0.091 





2 

53735.16000 

1.9 

0.087 





3 

53754.14737 

0.8 

0.051 


^to is the first time of dimming points that detected by CPA (in units of MJD). 

is the duration between the start time of eclipse and the end (in units of hours). 
is the depth of observed dimming points. 
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17.21 
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15.09 
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15.09 







53731.14064 


53742.42521 


53756.13038 


MJD (days) 

Fig. 11.— Zoomed-in light curve collection of 
eclipse candidates. Short vertical lines in the mid¬ 
dle of each panel mark epochs of eclipses. The 
numbers in the upper right denote the duration 
(At) between the start time of eclipse and the end. 


spectively), and has out-of-eclipse variability at 
the 2% level with period ~4.0631 days. After 
pre-whitening with the rotational modulation, our 
analysis indicates possible orbital period Porb = 
2.4277 days. 

V926 This object was known to be a long-period 
(16.5777 days) variable star. Our LS analysis 
supports the conclusion of previous study, but it 
has sinusoidal-shaped periodic signal in our light 
curve. The LS periodograin shows two power 




0.4 0.45 0.5 0.55 0.6 




Fig. 12.— Phased light curves of V985, V1483, 
V1731, and V1889. The mid-eclipse is set to the 
orbital phase of 0.5. The estimated orbital period 
is given in the upper right of each panel. 


peaks with periods similar to previous one (rank 
1 = 15.7291 days) and another shorter than one 
day (rank 2 = 0.9376 days). One sudden eclipse 
event was caught during its return to quiescence 
brightness. 


VI018 We caught two eclipses, but the star 
seems to vary all the time in a strange manner. 


5.1.4-- Grazing EB candidates 


By using the Plavchan algorithm (|Plavchan et al 
I 2 OO 8 I provided in NASA Exoplanet Archive Pe- 
riodogram Service]!, we identified two types of 
systems which can be confused with a transiting 
planetary system. One type is connected with 
the grazing EBs, and it shows V- or U-shape 
eclipses. The other type is the eclipsing binary 
systems with a large primary star, in which case 
the eclipse produces a flat bottomed light curve 
with small deeming. Figure [T^] shows the phased 
light curves of actual cases of these contaminants. 
The observed shapes of the eclipsing events resem- 


®NASA Exoplanet Archive Peri- 

odogram Service is available at 

http://exoplanetarchlve.ipac.caltech.edu/egi-bin/Periodogram/nph-simple: 
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ble transit-like signals, but depths of mid-eclipse 
are more than a few % of the total light (V985: 
5%; V1483: 14%; V1731: 24%; V1889: 10%). 


5.2. Multiperiodic variable stars 


Multi-periodicity is a common feature found in 
nearly all type of pulsating stars such as d Scuti 
stars, RR Lyraes, and Cepheids. Our catalog 
contains many examples of multiperiodic variable 
stars which deserve in-depth analysis. 

We performed multiple-frequency analysis of 
92 short-p eriodic variables with the program PE- 
RIOD04 ( Lenz fc Breged 2005h . looking for the 
presence of multi-periodicity. With an iterative 
pre-whitening procedure, discrete Fourier trans¬ 
formations are calculated until there is no mean¬ 
ingful change in the fit residuals. Each light curve 
was fitted using an equation of the form: 


y{t) =Ao + y] Aj cos(27r/j t - </)j) (9) 


1=1 


where Aq is a mean level and fj, Aj, and tpj are 
the frequency, amplitude, and phase for each suc¬ 
cessive peak found in the amplitude spectrum. We 
adopt a conservative approach in selecting the sta¬ 
tistical significant peaks from the amplitude spec¬ 
trum. In general, a signal/noise (S/N) amplitude 
ratio of 4.0 is a good criterion for independent fre- 
que ncies, equivalent to 99.9% certainty of var iabil- 
ity ( Breger et al ] l 199.31: iKnschnief et al ] ll997ll . For 


combination modes with known values, the S/N 
amplitude ratio of ^3.5 is used which is equivalent 
to 90% certainty. After removal of all significant 
peaks in the amplitude spectrum, the noise level 
was estimated on the residual amplitude spectrum 
over 5 c/d boxes around each frequency. The re¬ 
sulting light curve of one example multiperiodic 
variable is shown in Figure 1131 The multi-period 
model is in excellent agreement with the observed 
light curve. 

Table IASI summarizes all candidate frequencies 
identified in 92 multiperiodic light curves with 
their fitted Fourier parameters. All the detected 
frequencies are classified as independent modes 
(/i, / 2 , /a, ■ • ■): combination modes {afi ± bfj), 
and harmonic mode {hfj). 
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Fig. 13.— Example light curve of multiperiodic 
variable V1898. The fit of the 16-frequency solu¬ 
tion presented in the Table IA2I is shown as a solid 
curve. The full amplitude of this variation is about 
0.04 mag in r-band light curve. 


5.3. Aperiodic variability 

By using the method described in Section 3.1, 
we detected long-term aperiodic variability of 132 
variable candidates with time-scales similar to the 
total observing time span (^31 days). We visually 
inspected each light curve to confirm its variable 
nature, and this result is illustrated in terms of 
time sequence (Figure [H]). Due to the lack of pe¬ 
riodicity, the light curves of aperiodic variables are 
simply characterized by the level of photometric 
variability (Am), the dispersion about the mean 
level (iTm), and the reduced chi-square assuming a 
constant median value for the light curve (xD- 


1 (xi - x)'^ 

N-1^ 

i * 


( 10 ) 


16 


where at is the photometric uncertainty of each 
data point Xi, and N is the total number of mea- 
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Fig. 14.— Examples of aperiodic light curves 
showing variations about the mean, in addition 
to a relatively long-term variability on timescales 
similar to the total observing time span (~31 
days). 



('•-Oo 


H-K, 


surements. The full parameters are described in 
Table Edl These variables show relatively low lev¬ 
els of variability with Am < 0.03 mag (60%) and 
Am < 0.05 mag (80%), respectively. 

We used the combined SDSS and 2MASS 
photometry to identify basic properties of these 
sources. Since color ranges are much smaller in 
the near-infrared color-color diagram, the stellar 
locus is not as well dehned as in optical color-color 
diagram (Figure fT^ . We infer spectral types from 
Covey et al.l ([20071), in which stars in this sample 
cover spectral types from approximately A5 to M5 
and can be either main-sequence and giant stars. 
The CMD indicates that many of these aperiodic 
stars are held stars not associated with M37. 


The origin of aperiodic nature is difficult to 
explain. Possible speculations include that (i) 
some of these may be cases where long-period 
held giants are detected as low-amplitude photo¬ 
metric \ 2 JmWes_onEnifiscaffis_o|_weeks_ 0 £longer 
(e.g., Henry et al. 20001 : Ciardi et al.l 2011 1. (ii) 
some of them may be semi-regular variables ex- 
hibiting both fast and slow variabilities (e.g.. 


Androno v fc Chinaroval 1 201311. a nd (iii) as men¬ 


tioned in Hartman et al.l ( 2008bll . the rapid am¬ 
plitude variation may be dickering, nova-like out¬ 
bursts seen in many symbiotic variables (e.g.. 


Fig. 15.— Color-magnitude diagram and color- 
color diagrams of aperiodic variable candidates 
(black points). The color index distributions 
follow the l ocus of normal stars taken from 
ICovev et al.l ( 2007 1. The two lines show the lo¬ 
cation of the median stellar locus of dwarf (V) 
and giant (HI) stars having solar-metahicity. The 
cluster’s extinction is taken into account for SDSS 
color-color diagram. In order to derive AxjAy 
values for SDSS photometry system , we used the 
extinc tion coefficients compiled in Girardi et al.l 
(l2004ll . 


Gromadzki et al.l l2009l 2013 ). Further observa¬ 


tions of much longer duration will be necessary to 
properly characterize these aperiodic variables. 


5.4. Flare stars 

Earlier analysis of M37 data (Hartman et al. 
2008a) could not address hare events properly be¬ 
cause many data points were removed as outliers. 
Our new pipeline, with nearly 100% data utiliza¬ 
tion, offers opportunity to investigate hare stars 
more thoroughly. 

By CPA analysis, we found that 436 stars (625 
hare events) show haring activity characterized by 
a rapid hux increase and a slower decay. As ex- 
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Fig. 17.— Comparison between our and archival 
light curves of the two flare events. As shown in 
Paper I, new analysis provides more improved r- 
band data in terms of both photometric accuracy 
and data recovery rate. 


Fig. 16.— Color-magnitude and color-color dia¬ 
grams of flaring stars, detected in archival (left 
panels) and our new data (right panels), respec¬ 
tively. The red lines represent the m ain-sequence 


stella r locus with solar metallicity (jCovey et al 
I2nn7h . Most of these stars can be regarded as 
late-type main sequence stars (K4-M6) in the M37 
open cluster. 


pected, the location of these flaring stars on color- 
magnitude and color-color diagrams suggest that 
most of them are late-type main-sequence stars 
(K4-M6) in M37 open cluster (Figure [16]). Flares 
are also observed from earlier spectral types (AO- 
K3), but these are only a few and probably belong 
to field stars. 


The underlying mechanism of flares in stars 
cooler than about F5 is essentially the same as 
in the Sun. It is believed that these cool stars can 
generate a magnetic field through a dynamos at 
the interface betwee n the radiative core and the 
convective envelope ( Benz fc Giidell 201(TI 1. Thus, 
flaring is an inevitable result of strong magnetic 
field. But flares in hot stars, especially in A-type 
stars, are not easy to explain in the context of 


dynamo process because these stars should not 
generate strong magnetic fields in their thin con¬ 
vection layer and are incapable of storing larger 
energy in that region (lBalonall2012r) . 

The observed flare signature of one A-type star 
(V72) might be the result of binarity. If an A- 
type main sequence star has an active M dwarf 
secondary star, such binary star would look like a 
flare-generating A-type star. 


The number of flare events is increased by a fac¬ 
tor of two in the new analysis because our new cal¬ 
ibration procedure provides more accurate r-band 
light curves fFigure ITTI) . In Table fATl and Figure 
da we list the parameters of flare events discovered 
by our CPA analysis. They have a broad range of 
peak amplitude (Ampeak = 0.006 ^ 3.719 mag) 
and typical durations of a few minutes to several 
hours, irrespective of their spectral type. More 
detailed analysis of flare statistics and its relation 
to magnetic activity and stellar rotation will be 
presented in subsequent papers. 
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Fig. 18.— Various light curves of discovered flare events selected by color cut {{g — r)o < 1). Most of them 
show a sharp rise and subsequent exponential decay. 
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6. Conclusion 

We cataloged and investigated a variety of 
photometric variables in the field of M3 7 using 
our new light curves, which were obtained from 
archival imaging data of one-month long exo¬ 
planet transit survey with the 6.5-m MMT tele¬ 
scope. This new dataset allows a rare opportunity 
to explore different types of variability of short 
(^minutes) and long (^one-month) time-scales 
simultaneously. 

We find 61 eclipsing binary systems (EB), 92 
multiperiodic variable stars (mp), 132 aperiodic 
variables (A), and 436 flare stars (F), as well as 
several hundreds of rotating variables (R). Our 
work has increased the number of known vari¬ 
able stars in this field by 60%. The new vari¬ 
able catalog has been made available on the web 
(http://stardb.yonsei.ac.kr/). 
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Table A1 

New catalog of 2306 variable stars in M37 field 


VarlD 

StarlD 

(J2000) 

Period®' 

(days) 

Amp® 

(mag) 

(mag) 

r” 

(mag) 

g - r^ 
(mag) 

r — 
(mag) 

(mag) 

yd 

(mag) 

s - yd 

(mag) 

r 

(mag) 

J - 
(mag) 

H - Ks‘ 
(mag) 

VarType^ 

VI 

055220.39-^323319.65 




13.677 

0.285 

0.191 


13.406 

0.609 

12.447 

0.270 

0.063 

EB 

V2 

055216.58-^322815.32 



16.792 

14.840 

0.564 

0.290 

14.471 

14.653 

0.688 

13.539 

0.287 

0.134 

EB 

V3 

055233.01-^323241.91 

0.4225 

0.320 

17.982 

15.789 

0.560 

0.329 

15.284 

16.051 

0.728 

14.267 

0.355 

0.012 

EB 

V4 

055253.26+323301.47 

0.5582 

0.235 

18.079 

15.634 

0.762 

0.477 

15.093 

15.948 

0.938 

13.951 

0.497 

0.101 

EB 

V5 

055300.63+322450.81 

0.2788 

0.310 

18.036 

15.915 

0.459 

0.481 

15.235 

16.184 

0.585 

14.421 

0.359 

0.084 

mp 

V6 

055150.53+323234.70 

0.1098 

0.343 

18.218 

15.973 

0.591 

0.173 

15.622 



14.739 

0.342 

-0.009 

mp 

V7 

055239.09+323631.27 

0.3577 

0.433 

19.858 

17.252 

1.039 

0.056 

17.065 

17.890 

0.815 

15.779 

0.446 

0.196 

EB 

V8 

055234.32-^323218.70 

0.1195 












u 

V9 

055214.91-^322441.35 

0.0909 


14.873 

13.233 

0.137 

0.128 

13.114 

13.298 

0.370 

12.348 

0.171 

0.042 

p 

VIO 

055200.49+323648.19 

0.9432 

0.075 

16.844 

14.151 

0.928 

0.436 

13.493 

14.762 

0.531 

12.356 

0.567 

0.131 

p 


Note.— Table IaTI is published in its entirety in the electronic edition of Astronomical Journal. A portion is shown here for guidance regarding its form and content. 

^'Dominant period and its full amplitude. All blank fields are non-periodic variables. In cases where the perio dicity analysis is affecte d by a large scatter of light curves due to 
photometric bias or source blending, we used the already known period in that source from the previous catalog iHartman et ani2008tJ) . 

and z magnitudes from the Sloan Digital Sky Survey (SDSS) DR7 dAbazaiian et ani200^ . 

^g, r, and i magnitudes from [Hartman et al.l d2008al) and SDSS DR7 catalog, 
and V photometry from fKalirai et al.| d200lD . 

® J, H, and Kg photometry from the Two Micron All Sky Survey (2MASS) Point Source Catalog dCutri et al.|[^03[) . 

^Type of variability. The following types of variability were assigned: F = flare stars, P = pulsating variables, mp = multiperiodic variables, R = rotating variables, EB = 
eclipsing binary systems, A = aperiodic variables, and : (or var:) = variable candidates. We also use the two cases: nonvar = non-variables listed in Table [Tl U = These variable 
objects fell on chip gaps and were not observed. 




















Table A2 

Candidate frequencies identified in 92 multiperiodic light curves 


No. 

VarlD 

Star 

(J2000) 

r 

(mag) 

Frequency 

(cd-l) 

Frequency 

name 

Ar^ 

(mmag) 

<t> 

(deg) 

S/Nb 

Note‘s 

I 

V5 

055616.79+322514.20 

15.915 

3.5872 

fi 

310.1 

0.3331 

68.9 






7.1768 

f2(=2fl) 

35.4 

0.2906 

11.3 






10.7619 

f3(=3fl) 

22.8 

0.0418 

10.3 






5.6768 

h 

19.5 

0.2369 

5.7 






14.3473 

fs 

13.8 

0.9990 

7.1 


2 

V6 

055506.89+323303.20 

15.973 

9.1040 

fl(=f3-f2) 

343.4 

0.5433 

172.7 






18.2081 

f2(=f3-fl) 

89.9 

0.2201 

53.6 






27.3122 

f3(=fl+f2) 

27.6 

0.9237 

19.7 






8.7545 

h 

15.4 

0.5801 

8.1 






11.9319 

h 

14.5 

0.5557 

6.3 






36.4178 

f6(=2f2) 

10.1 

0.0527 

10.4 






14.9868 

f? 

00 

bo 

0.9887 

4.5 






16.8695 

fs 

7.9 

0.7280 

4.5 


3 

Vll 

055519.81+323541.23 

15.245 

8.4261 

fl 

79.0 

0.9052 

84.8 






16.8561 

f2(=2fl) 

7.6 

0.6488 

8.3 






17.7044 

f3 

8.0 

0.5895 

9.1 






16.4751 

f4 

6.9 

0.5617 

7.3 






8.8787 

fs 

4.9 

0.4440 

5.2 






12.0801 

fe 

4.8 

0.8357 

5.0 






9.2640 

f7 

5.9 

0.6453 

6.4 






10.6455 

fs 

5.4 

0.0597 

6.0 






13.5604 

fg 

4.8 

0.8005 

4.9 






16.2578 

fio 

5.1 

0.6379 

5.3 






8.0470 

fii 

4.5 

0.8268 

4.7 






2.3825 

fl2 

5.2 

0.7274 

3.6 






4.3812 

fl3 

4.3 

0.3942 

4.2 






14.9483 

fl4 

3.9 

0.8404 

4.0 


4 

V12 

055600.99+323041.25 

15.873 

12.7176 

fl 

14.4 

0.0994 

8.1 






7.8896 

f2 

13.8 

0.4901 

6.0 






12.7698 

f3 

10.5 

0.0175 

6.0 

close to fl 





14.4564 

f4 

9.6 

0.0884 

6.5 


5 

V13 

055527.32+322543.01 

16.511 

10.0341 

fl(=fl2-f2) 

23.0 

0.4259 

137.9 






10.4706 

f2(=fl2-fl) 

13.0 

0.8171 

77.9 






17.1389 

f3 

8.7 

0.9050 

51.9 






14.0376 

f4 

6.1 

0.6801 

33.6 






13.0968 

fs 

3.7 

0.1837 

20.9 






12.8019 

fe 

3.1 

0.2056 

17.4 






15.0841 

f7 

2.1 

0.7887 

11.6 






13.6263 

fs 

1.3 

0.1260 

7.1 






11.1091 

fg 

1.4 

0.2454 

8.7 






27.1718 

fio 

1.0 

0.7296 

6.7 



Note.— Table is published in its entirety in the electronic edition of Astronomical Journal. A portion is shown here for 

guidance regarding its form and content. 

^We doubled amplitudes derived by PERIOD04 since this program gives half of the full amplitudes. 

'^The signal-to-noise ratio (S/N) was calculated by PERIOD04. Each noise level was computed from the prewhitened periodogram 
as a running mean over boxes of 5 cd“^ in frequency. 

"^We note the frequency values at which the close pairs occur with a separation less than 0.1 cd“^ 
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Table A3 

Long-term 132 aperiodic variable candidates with timescales of >1 month 


No. 

VarlD 

Star 

(J2000) 

r 

(mag) 

(mag) 

Arn° 

(mag) 


Nlc 

Note‘s 

1 

V58 

055122.14+322942.40 

16.742 

0.009 

0.096 

443.5 

4730 


2 

V70 

055123.93+324357.22 

15.837 

0.003 

0.017 

14.0 

4706 


3 

V75 

055124.73+324312.36 

18.360 

0.004 

0.028 

1.8 

4730 

13.260 

4 

V87 

055125.42+323600.45 

17.564 

0.002 

0.017 

1.6 

4730 

21.199 

5 

V107 

055127.04+323225.57 

15.458 

0.002 

0.010 

7.7 

4592 


6 

V1I4 

055128.19+322931.14 

17.854 

0.004 

0.029 

3.4 

4730 


7 

V131 

055129.26+323235.44 

18.515 

0.005 

0.029 

1.4 

4730 


8 

V147 

055130.48+324053.05 

18.348 

0.005 

0.033 

1.9 

4730 

7.709 

9 

V157 

055131.24+323059.61 

17.688 

0.008 

0.046 

12.9 

4730 


10 

V163 

055131.87+324233.83 

17.482 

0.004 

0.026 

5.1 

4730 



Note.— Table is published in its entirety in the electronic edition of Astronomical Journal. A 

portion is shown here for guidance regarding its form and content. 

^cTm is the measured standard deviation of Scr-clipped light curve. 

^Am is the peak-to-trough amplitude of Scr-clipped light curve (Am = mmax — '^min) 

^The reduced chi-square ixt) of the magnitudes is to test the significance of variability for individual 
objects against the null hypothesis of no variation. 

"^For the objects that were pr eviously recognized as periodic variables, we fill the column with the 
period value using the catalog of lHartman et ^ ll2008br) . 
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Table A4 
436 FLARE STARS 


No. 

VarlD 

Star 

(J2000) 

r 

(mag) 

FlarCevent 

(No.) 

j. a 

^peak 

(days) 

^rUpeak 

(mag) 

rO. 2 '' 

(hours) 

Note'* 

1 

V42 

055120.47+322200.75 

20.690 

1 

53742.40307 

0.093 

0.17 






2 

53742.46247 

0.135 

0.08 


2 

V50 

055121.46+322856.90 

19.990 

1 

53730.34740 

0.221 

0.55 






2 

53733.19093 

0.176 

0.21 






3 

53735.19767 

0.039 

0.49 






4 

53737.33376 

0.188 

0.30 


3 

V69 

055123.45+322948.83 

20.809 

1 

53732.22482 

0.061 

0.32 






2 

53732.35513 

0.133 

1.11 






3 

53742.35427 

0.281 

0.23 


4 

V72 

055124.44+324407.71 

15.143 

1 

53730.32306 

0.010 

1.47 


5 

V76 

055124.78+323440.56 

18.368 

1 

53728.24480 

0.026 

0.17 






2 

53742.44376 

0.117 

0.24 






3 

53753.25302 

0.139 

0.16 


6 

V78 

055124.82+322937.32 

20.017 

1 

53726.39533 

0.046 

0.10 

Var 

7 

V86 

055125.41+323800.03 

18.086 

1 

53726.27563 

0.041 

0.58 

Var 





2 

53727.30272 

0.025 

0.46 


8 

V91 

055125.71+323106.38 

20.553 

1 

53734.23990 

0.072 

0.20 

Var 





2 

53756.08268 

0.079 

0.60 


9 

V104 

055126.78+324024.73 

17.971 

1 

53726.43908 

0.030 

0.29 


10 

V112 

055127.86+322854.79 

22.181 

1 

53727.19210 

0.360 

0.58 

Var 


Note.— Table lA4l is published in its entirety in the electronic edition of Astronomical Journal. A portion is shown 
here for guidance regarding its form and content. 

^tpeak is the time at flare peak. 

^Arripeak is incremental magnitude at flare peak. 

‘^We define the flare duration (tq, 2 =^ 0 .2,rise + ^o. 2 ,decay) ^s the time-scales of flare events at 80% of its peak level 
(Paper III). 

^These stars show periodic brightness variations with a near or distorted sinusoidal shape (Var). 
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